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Abstract 

In this paper, we look at the effect of volatility clustering on the risk indifference 
price of options described by Sircar and Sturm in the paper ‘From smile asymptotics 
to market risk measures’ [l2]. The indifference price in [12] is obtained by using 
dynamic convex risk measures given by backward stochastic differential equations (BS- 
DEs). Volatility clustering is modeled by a fast mean-reverting volatility in a stochastic 
volatility model for stock price. Asymptotics of the indifference price of options and 
their corresponding implied volatility are obtained in this paper, as the mean-reversion 
time approaches zero. Correction terms to the asymptotic option price and implied 
volatility are also obtained. 

Keywords: risk measures, indifference price, implied volatility, volatility clustering, 
comparison principle. 


1 Introduction 

In an incomplete market, there are several ways of pricing options. The arbitrage-free method 
of option pricing, where the option price is given by the expected value of the discounted 
payoff under a risk-neutral equivalent martingale measure, is widely known and studied (see 
for example [1]). Another method is the indifference pricing method. Indifference pricing of 
options using utility functions have been studied extensively in 0 . 0 . 1111 . etc. Later it was 
observed that the same idea of indifference pricing could be extended from utility functions 
to dynamic convex risk measures, see [2]. Sircar and Sturm in na derived a nonlinear 
partial differential equation (PDE) characterizing the indifference price of put options given 
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by dynamic convex risk measnres. In their paper, they used the residual risk measure after 
hedging for pricing options; the residual risk measure was given in terms of a BSDE. They 
also obtained the implied volatility, corresponding to this indifference option price, as the 
viscosity solution of a nonlinear PDE. 

The significance of the results in [12] is that, via the indifference pricing scheme, the 
market risk, reflected in the implied volatility skew, can be related to convex risk measure 
theory. Typically, convex risk measures are defined abstractly via BSDEs (see 0 )- So the 
right driver in the BSDE which gives a good risk measure is hard to determine. Using the 
indifference pricing of Sircar and Sturm, we can calibrate the driver from the market implied 
volatility data. In their paper, na, small-maturity asymptotics of the implied volatility 
yielded simple formulas. However, in general, a closed form solution to the non linear PDE 
in [12] is hard, if not impossible, to hnd. In this paper, an effort is made to give some 
meaningful, simplihed formulas for the implied volatility surface, by considering the effect of 
volatility clustering. 

It is believed that market volatility fluctuates frequently between high and low periods. 
While volatility cannot be observed directly, this “clustering” behavior is estimated from 
observed stock prices. To model stock prices subject to this clustering behavior of volatility 
we use stochastic volatility models where the volatility is a fast mean-reverting ergodic 
process. Such stochastic volatility models with fast mean-reverting volatility were found to 
be a good £t for stock price data (see chapter 4 of m)- The question of interest in this 
paper is, how does fast mean-reversion in volatility affect option prices? As the rate of 
mean-reversion increases, the long-run behavior of the ergodic volatility process manifests. 
Consequently, the effect of volatility gets averaged with respect to the invariant distribution 
of the ergodic volatility process. The analysis of fast mean-reversion of volatility on option 
prices has been studied in the case of no-arbitrage pricing, see chapter 5 of [1] . In this paper, 
we will look at the effect of fast mean-reversion of volatility on the indifference prices of 
options given in na. 

In the no-arbitrage pricing case, option prices are given as the solution of linear PDEs. For 
small mean-reversion time, denoted by the parameter e, this leads to a singular perturbation 
problem. In |3] , assuming the volatility process is an Ornstein-Uhlenbeck process, the option 
price is expanded in powers of e and the asymptotic option price, when e —)■ 0, is obtained 
together with a correction term of order ^/e. The corresponding corrected implied volatility 
is also obtained. 

In this paper, unlike the no-arbitrage case, the option price is not the solution of a linear 
PDE. In fact, this indifference price is given in terms of the solutions to BSDEs. Initially, 
finding the asymptotic indifference price appears to involve averaging of BSDEs. However, 
there is a difficulty: Sircar and Sturm use quadratic drivers in their BSDEs which are not 
Lipschitz and hence we cannot use the established stability results for BSDEs (for example 
in m)- We instead avail of the nonlinear Feynman-Kac formula derived by Pardoux and 
Peng in uni, and express the option price in terms of solutions to nonlinear PDEs; thus 
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converting this into a problem of averaging nonlinear PDEs. 


The formal derivation of the corrected option price and implied volatility, in sections 3.1 


and 1^ follows the same line of reasoning seen in [3]. The difference in method lies in the 
rigorons proof of accnracy of the corrected asymptotic formula, which is in section |3.2[ For 


the rigorous proof, we use the maximum principle to bound the solution of the nonlinear 
PDF, which gives the option price, within o{^/e) distance of the corrected asymptotic option 
price formula. 

The purpose of this paper is two-fold: on one hand, we look at the effect of an observed 
phenomenon, viz., volatility clustering, on option pricing; secondly, the averaging of the 
volatility leads to a simpler formula for implied volatility, which makes calibration easier. The 
dependence of the corrected implied volatility formula on the risk parameters can potentially 
be used to calibrate risk measures from the market implied volatility data. 

It should be interesting to see if this work extends to other risk measures besides those 
given by BSDEs. This however falls outside the scope of this paper and will be considered for 
future work. In general, regardless of risk measure used, if the security price can be expressed 
as a solution of a PDF, we can potentially use the method in this paper for averaging out 
the effect of fast mean-reverting volatility. The paper is organized as follows. In section 
we recall the indifference pricing of options from [12] and introduce the nonlinear PDF 
that gives the indifference price of put options. Section 3.1 has heuristic calculations for the 
corrected asymptotic option price. The main result and the rigorous proof of accuracy of the 
obtained corrected asymptotic option price formula is in section 3^ The corrected implied 
volatility formula is obtained in section 


2 Preliminaries 

We begin by introducing the stochastic volatility model for stock price. 

2.1 Stochastic volatility model. 

Let denote stock price at time t, where the parameter e refers to the mean-reversion 
time scale for volatility. As e approaches 0, the speed of mean-reversion, 1/e, increases. 
Let (fl, P) denote the probability space on which Se,t satishes the following stochastic 
volatility model. 

=b{Y,^t)S,,tdt + a,{Y,^t)S,,tdW,^^\ 0 < f < T, 

dY.^t + y/l 0 < f < T, 

e ye 
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(la) 

(lb) 

(l c) 









where |p| < 1, and are independent Brownian motions on (h2, P). We make the 

same assnmptions on the stochastic volatility model as in |T2] . For the reader’s convenience 
we recall these assnmptions. 

Assumption 2.1. We assume that 

1. ai,(T2 e where C;;/{R) is the space of differentiable functions with locally 

Holder-continuous derivatives with Holder-exponent (3 > 0. 


2. Both (Ti and <72 are bounded and bounded away from zero: 

0 < Cl < (Ji < Cl < cxD, and 0 < C 2 < <72 < C 2 < oo, 


3. b ^ and b is bounded. 

Let B denote the infinitesimal generator of the Y process when e = 1. Then, for / G 
C2(M), 

Bf{y) := (m - y)dyf{y) + ]^(yl{y)dlyf{y). (2) 

By the general theory of 1-D diffnsions (see Karlin and Taylor [8], page 221) it is easy to see 
that there exists a nniqne probability measnre 


^{dy) = Z 


-1 


exp 


y 2{m-z) ,^ 
0 


(^2iy) 


dy, 


(3) 


snch that f Bf{y)7i{dy) = 0 for all / G Z is the normalizing constant, so that 

J 7i{dy) = 1. The invariant distribntion of Y given by ([^ plays an important role in the 
following analysis. 


2.2 Indifference option price. 

We consider a European put option with maturity time T and strike price K. The indifference 
price of this European put option at time t, Pe{t, x, y), is given in [12] in terms of risk measures 
as follows. 

P,{t,x,y) = R^^t-Re,t, (4) 

where {R, Z.) and {R., Z.) are respectively solutions of the following BSDEs 

R = - fW,\ ziWs - z'fyw-f >; (5) 
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and 


R,f = -(A' - S,,3.)+ - /(Zi’i, Zl?)* - Z«<iW'i« - ZiJdW'f >. (6) 

The function / in the above BSDEs satishes the criteria for admissible drivers (see Dehnition 
2.3 in ini) to ensure the solvability of the BSDEs. In this paper, we will consider a specihc 
family of admissible drivers called distorted entropic risk measures which were introduced in 
[12] (see section 3.1 of [I2|). This class of drivers has the following form: 

9'^’^{zi,Z2) := I {{zi + riZ 2 y + zj) , 

and is parametrized by two parameters: the risk aversion parameter 7 > 0 and the volatility 
risk premium rj. When p = 0 the driver reduces to the classical entropic risk measure 
whose level curves are circles with radius depending on the risk aversion parameter, 7 . By 
introducing p, we distort this circle into an ellipse (if \ri\ < 1). Under hedging, the risk 
measure gets adjusted and is now given by a BSDE where the driver g'^’Z is transformed to 


g^’'^{zi, Z 2 ) = inf (g^’^(zi + ai(|/)z/, Z 2 ) + b(y)iy) 

b{y) b^jy) 7 2 Vb{y)z 2 

^Vi(?/) 2'yal{y) 2^^ ai{y) 


Henceforth we will take this to be our driver, i.e. dehne 


f{Zi,Z2) 


b{y) b^jy) 
^Vi(?/) 2'jaf{y) 


7 2 , Vb{y)z 2 
2 "^ ^ o,{y) ■ 


It is important to note that this driver is not Lipschitz, as it is quadratic in Z 2 

Instead of stock price S^^t, we will work with the logarithm of stock price normalized by 

the strike price: ■= In (%^ )• With this change of variable, ([^ becomes 


= (^6(y;,i) - dt + 0<t<T, (7a) 

dY,t + y/l - p^dW^^'^), 0<t<T, (7b) 

e ^/e 

(W,o,b;,o) = (x, 2 /). (7c) 


Sircar and Sturm in na use the generalized Feynman-Kac formula, given by Pardoux 
and Peng in [10], to describe the solutions of the forward-backward SDEs, and R^^t, in 
terms of solutions of the nonlinear PDE ([^ below. Before introducing this PDE, we will 
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first make another change of variable t -.= T — t, which gives the time to maturity. Define 
the differential operator by 


^ 2 / \ o2 ^ 2 / \ o2 P 


^e9 ■=i^(^i{y)d^x9 + w:^2{y)dyy9 + —(^i{y)(^2iy)d^y9 


2 e 


'xy^ 


1 . 

-[m - 

h\y) 

2-ial{y) 


y) 


JLm^\ 

a,{y) ) 2 


( 8 ) 


+ 


7(1 „^2 9\/^- p^Ky)(^2{y)dy9 


2 e 


-<yi{y){dy9y- 






for g G Let and denote solutions to the PDE 


drU = C^u, (9) 

with the initial conditions Ue{0,x,y) = 0 and Ue{0,x,y) = —[K — respectively. By 

Theorem 2.9 in [12], we get the put option price 

Pe{T,x,y) = Ue{T,y) - u^{T,x,y). ( 10 ) 


As mentioned in na, by Ladyshenskaya et ah jS] Theorem V.8.1], and are unique 
bounded classical solutions to the semilinear parabolic equation ([^, with bounded derivatives 
in [0, T] X M X M. 

Remark 2.1. Since the coefficients of the PDE in 0 and the initial condition of u are 
x-independent, we get u to be x-independent. 

Remark 2.2. While we only consider European put options, the results in the paper extend 
to any other option with bounded and continuous payoff. 


3 Asymptotic option price 

We will begin with heuristic arguments for obtaining the asymptotic option price and cor¬ 
rection terms to the asymptotic price. The heuristic calculations follow along the same lines 
as the no-arbitrage option pricing case, seen in | 1 | 


3.1 Heuristics 

Assume the following expansions of and P^ in powers of -y/e: 


We = Mo + a/cmi -h eu2 + -1- ..., 

(lla) 

u^ = Uo + y/eui + eu 2 + -1- ..., 

(11b) 

Pe = Pq + + eP2 + . . . . 

(11c) 
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We define the following differential operators: for g E xM.), 

^o9{x, y) : = y) - \al{y)d^g{x, y) - 


b\y) 


2^afiy)' 


^i9ix,y): = pai{y)a2iy)dlyg{x,y) 


b{y)(^2{y) 

.A , 7(1 Kl{y){dyg{x,y)f . 


- (p + p\/l -p^) -^b^r7^9yg{x, y), 
^2^9{x,y)-. = Bg{x,y) + 


Note that the operator is nonlinear. Equation (|^ can be rewritten as 


drUe = CeUe = CqU^ H-p/ilMe H- 


( 12 ) 


3.1.1 Leading order term 


Using the expansion (11a) in (12) and collecting terms of order 1/e we get 

7(1 


Bun = 


-a': 


:{y){dyUo) , 


which is satisfied if Uq is y independent. We will thus assume Uq{t,x) is independent of y. 
Using this and collecting terms of order I/a/c we get 

Bui = 0 , 


which is satisfied by taking ui independent of y. As Uo{t,x) and Ui{t,x) are both y inde¬ 
pendent, terms of 0(1) in equation ([^ satisfy 

drUQ{T, x) = CqUo{t, x) + BU 2 . 

Thus, U 2 must satisfy the Poisson equation 


Bu 2 = drUo{T,x) - CqUo, 

which has a solution provided the following centering condition holds: 

„ 1 _ 2 < a 2 1 _ 2 / 1 b‘^{y) 

OrUo = -CTI d^^uo - -ai {y)d^uo - - - 

2 2 ^7c^i(P) 


(13) 


(14) 


Here ai^ and ^ 5 ^ denote the average of the terms crf(y) and respectively, with respect 


o-ilU ■ 


to the invariant distribution, tt, of the Y process. Similarly, uq satisfies equation (14). 
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The initial conditions for and give the initial conditions for uq and uq respectively, 
i.e. Mo( 0 ) = 0 (ho is independent of both x and y, so only a function of r) and Mo( 0 ,a;) = 
— [K—Ke^Y^. Observe that the hrst order approximation term to the option price, Pq{t, x) = 
Uq{t) — Uo{t,x), satisfies the equation 

dtPo = ^^^dl^Po - ^P^^{y)dxPo, PoiO, x) = [K - Ke^'^, 

which is simply the equation for the Black-Scholes put option price, PBs{.'r,x;ai), with 
volatility parameter dq. Therefore, 

Po{r,x) = PBs{r,x-,a^) = KN{-d 2 ) - Ke^N{-di), (15) 

where N{z) = di = and ^2 = ^ ■ Observe that 

uo{t,x) = Po{t,x) - (16) 

27 


3.1.2 Higher order terms 

Define jCq := Cq + — ^aldx- So Cq is a linear operator. Equating terms of 

order ^/e in (|^ gives the equation 

drUi = CqUi + C1U2 + BU'i. 

Thus M 3 is the solution of the Poisson equation 

Bu^ = drUi — CqUi — CiU2- (17) 

provided the following centering condition is satished: 

drill = CqUi + ClU2- 


Henceforth, a line over any term indicates averaging of the term with respect to the invariant 
distribution of the Y process given in (|^. Since ui is independent of ?/, the above equation 
becomes ^ ^ 

drUl = -oY^dlr-Ui - -oY^dxUi + £iM 2 - (18) 

To simplify the right hand side of the above, recall that U 2 is the solution of the Poisson 
equation in (13). Together with (14), we see that U 2 is the solution of the following equation, 

b^{y) 


Bu2{T,x,y) = 


27 fM ■ 2^ + ^ Kto) --iq (d.u„{T,x)-alMrx)) 














Let 01 (y) and 02 (y) denote the solutions of 


BMy) = 




27af(2/) 27af(?/) 


and 


respectively. Define 


^02 (l/) = 2 K(l/) -^ 1 ^) 


Then 


t/2(r, x, I/) := 0i(?/) + 02(2/) (5^Mo(r, x) - dl^uo{T, x)) . 
M2 (d X, y) = f/2(r, X, y) + F(r, x) 


where F{t,x) will be determined later, see (27). 

Remark 3.1. Note that the function U 2 {T,x,y) satisfies equations (13) and 
We compute 


(19) 

( 20 ) 


C 1 U 2 


= p(Ti{y)a2{y)(t)'2{y) [dl^uo{T,x) - dl^^u^ir^x)) 


- [p + hVl -P^) {dxUo{T, x) - dl^uo{T, x)) 


Substituting this in (18) we get 


drUl 


= - ]^(jl{y)d^ui - Adl^^uo{T,x) 

+ (A + B)dl^ufiT, x) - BdxUo{T, x) - A, 


where 


and 


A = pai(2/)cr2(2/)02(2/), 


A 


(p + h\/l -P^) 


h{y)a2{y) 

^liy) 


01 (p) 


B 


(p + 7V^l -p2) 


h{y)(T2{y) 


02 (p)- 


( 21 ) 

( 22 ) 

(23) 
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It is easy to verify that the solution to the above equation is 

Mi(r, x) =T -Ad^^^uoir, a:) + (A + B)d^^Uo(T, x) - Bd^Uoir, x) - A 
By a similar argument, we get 


Mi(r) = T -Adl^^uoir) + (A + B)dl^uo{T) - Bd^uoir) - A 


= —At. 


Thus, 


Pi(r, x) = T [-Adl^^Po{T, x) + (A + B)dl^Po{T, x) - Bd^Poir, x)] 


3.1.3 Terms of order e 

Using the asymptotic expansion of in ([^ and collecting terms of 0(e), we get 


drU2 = jCoU 2 + ClUs + BU4 + 
Thus U 4 is the solution of the Poisson equation 


7(1-P^) 


^l{dyU 2 f. 


Bu 4 = drU2 - toU2 - P1M3 “ ^ ^ ^ CFl{dyU2Y, 


provided the following centering condition holds: 


5t'W2 ~ BqU 2 — PlU'i 


7(1 -p 2 ) 


-a. 


l{dyU 2 y = 0 . 


The unknown function F{t,x) in (20) is determined by the above equation (27). 


3.2 Accuracy of corrected asymptotic price formula 


( 24 ) 


(25) 


(26) 


(27) 


Theorem 1. Let P^ denote the put option price given by (10), and let Pq and Pi he as 


defined in (15) and (25), respectively. Then, under assumptions 2.1 


\Pe{r,x,y) - {Po{t,x) + ^fePl{T,x))\ < O(-eloge). 

Proof. We use a maximum/comparison principle argument to prove this result. 
Define 

My) = Ay - + De, 

where Df > 0 is a constant chosen such that 


(28) 


My) 


> \u2iT,x,y) \ + \/e \u3{T,x,y)\ 


(29) 
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for all {T,x,y) G [0,T] x M x 
most logarithmic growth in y (see (16), (24), 
functions of r and x (see (19) and (16)). It is easy to check that = 0{~ loge), see Lemma 
A. 2 | in the Appendix. 


It is possible to hnd such a. D^, as U 2 and M 3 have at 
, (p)^ and Lemma A.l) and are bounded 


Choose a constant C > 0 large enough so that 

C A 2 ( 1 / — mY' > drU2 — -^ 0^2 ~ -^ 1 ^ 3 . (30) 

It is possible to choose such a C since U 2 and M 3 are bounded functions of r and x and have 
at most logarithmic growth in y (see (16), (24), 

Dehne 


.0 and Lemma [a. 1[ ). 


■= Uoij^x) + i/eMi(r,x) + eU2{T,x,y) + x,y) - ei/j^y) - e{C + cI)t, 

and 

M,(r, x, y) := Mo(r, x) + \/eMi(r, x) + eU2{r, x, y) + e^/^M3(r, x, y) + eip^v) + YC + 


where C 2 is the upper bound on ct 2 in Assumption 2.1 , Mo,Mi, U 2 and Ye are defined in 
and (28) respectively. 


We compute 

drU^ — =drUQ + y/kdrUi + edrU2 + e^^'^drU^ — e{C + cl) 
— CqUq — \feC,QU\ — e£of^2 ~ 

— \feC,\U 2 + \f^C-\Ye ~ 

7(1-P>?(1/) 


EU 2 - \^Bu^ + BYe - e- 


{dyU2 + y/edyU 3 - dyYe) , 


where we have used the independence of y in mq and Mi to eliminate some terms. Using 
Remark 3.1, (13), (17) and BYe < —e2(|/ — m)^ + ec^, we get 


< — e{C + 2 {y — m)^) 


+ e' 


3/2 


- 2 (p + hVl -P^) 


drU2 — C0U2 — C1U3 


Ky)(^ 2 {y) 

^liy) 
7(i-P>?(i/) 


{y-m) + drU 3 - C 0 U 3 

{dyU 2 + yfedyU 3 - e 2 {y - m))' 


The term —e{C + 2 {y — m)Y dominates the other terms on the right hand side of the above 
equation, uniformly, for small enough e. To see this, hrst look at the terms of 0(e^/^). The 
hrst term grows linearly in y as b and a 2 are bounded and ui is bounded away from 0. The 
terms that depend on M 3 grow at most logarithmically in y and are bounded in r and x. 
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This can be seen from (16), (@, (0, (@ and Lemma |A.l[ Thus the terms are 

dominated by e{C+2{y—m)‘^) uniformly for all r, x, y, for small enough e. We now turn to the 


remaining term of 0(e) in the third line of the above inequality. By choice of O in (30), this 
term is dominated by e(0 + 2(?/ —m)^). Therefore the term —e{C + 2{y — m)‘^) dominates the 
other terms for small enough e making the right hand side of the above inequality negative. 
Hence, 

drU^{T,x,y) - C,u^{T,x,y) <0, W{t, x,y) e {0,T) x R x R, (31) 

for small enough e. Similarly, it can be shown that 


drU^{T, X, y) - X, y) > 0, V(r, x, y) e (0, T) x M x 


(32) 


for small enough e. 

By the maximum principle, we can show that < u^:, as follows. The goal is to 

prove that at the point of maximum of we have < u^. For this we will first need 

to ensure the maximum of is attained and to this end we first perturb our function 

slightly. Define 

' 5 (yrT^+ ^ 


T-t 


where 0 < 5 1. Note that the first and second order derivatives of \/l + x^ are bounded. 


Therefore the strict inequalities (31) and (32) still hold if is replaced with for small 
enough 6 . Because of the term —'il)^{y) in which dominates the growth in ?/, —>■ —cx) 

when \y\ —)■ oo. Also observe that is bounded in r and x, so the perturbation of by 
—5 (a/ 1 + x^ + j^) ensures nf —)■ —oo when either |x| —)■ cxd or r —)■ T. Since is bounded, 
we see that ^ —oo when either \y\ —)■ cx), |x| —)■ cxd or r —>■ T. Therefore, nf — 

must attain its maximum at some hnite point, say (ro,xo,|/o), in the domain [0,T) x M x M. 
For (tq, Xq, yo) in the interior of (0, T) x M x M, we have 


0 = drU,(T„,XQ,y„) - £,u,(To,Xo,yo) < Sru‘(Ta,Xo, Po) - ■C,S<f(To,io. !/l)) < 0; 


(33) 


which is a contradiction. Therefore the point of maximum must occur at the boundary r = 0. 
Recall that at r = 0, Mo(0,x) = u^(0,x,y) = —[K — and Mi(0,x) = 0 from ( pi) ). 

Then, by the choice of in the dehnition of V’e he. (29), we get Mf(0,x, y) < Me(0,x, ?/), for 
all (x, I/) G M X M and thus uf < everywhere, for all 5 1. Taking the limit at 5 0, we 

have the desired comparison 


u^iT,x,y) < u^{T,x,y), V(r,x,|/) e [0,T] x M x 

By a similar argument, we get 

Ue > ne(r, x,?/), V(r,x,|/) G [0,T] x M x M. 


(34) 


(35) 
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Putting inequalities (34) and (35) together and using (29), we get 


W{t,x,y) - {uq{t,x) + y/eui{T,x))\ < + e{C + cI)t, 

for all (r, x, y) G [0, T] x M x M. Recall that the constant in the dehnition of V’e is of order 
— log(e) and that V’e has a term with quadratic growth in y, which gives us 

\u,{t,x,y) - iuo{T,x) + y/eui{T,x))\ < O(-eloge) + 0 {e‘^)y‘^. 

We can repeat the same argument for which leads to the desired result 

\P^{t,x,y) - {Po{t,x) + ^/ePl{T,x))\ < O(-eloge) + 

for all {T,x,y) G [0,T] x M x M. □ 

4 Implied Volatility 

In theory, we could use the corrected option price formula for calibration. However, option 
prices are typically quoted in terms of their implied volatility and so it is useful to derive a 
corrected implied volatility formula corresponding to the corrected option price. Let denote 
the implied volatility corresponding to the put option price P^. Recall that x = ln{S/K) 
and the Black-Scholes formula for put option prices with volatility a is given by the formula 

Pbs{p x; a) = KN{-d2) - iPe"iV(-di), 

where N{z) = ^ and 4 = 

By dehnition, the implied volatility Jg is obtained by setting 


PBs(T,X',I,(T,x,y)) = P,(T,X,y). 


(36) 


From (15), we know that the leading order term of the put option price Po{t, x) = Pbs{x, x; oF). 
Thus, the leading order term of the implied volatility is simply cfi- To obtain a correction 
term for the asymptotic implied volatility we will expand R about dq in powers of -y/e as 
follows 

4('r, x,y) = d\ + A/e/i(r, x, y) + ehij, x, ?/) H-. (37) 


Using the expansion of U in (36), we get 


FB5(r, x;(Ti) + A/e/i^^^(r,x;(Ti) H- 

= ^o(u x) + yePi(r, X, y) H-. 


(38) 
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Thus the correction term to the implied volatility is 


h = Pi{T,x,y) 


QPbs / _ 


n -1 


Differentiating the Black-Scholes formula, we get 

dPBs 

da 


Substituting the formula for Pi in (39), we can rewrite (37) as 


h = 0-1 - ve 


Ki/re '^ 2/2 
Bd:^Po{T,x))] + o(^/e), 


x) + {A + B)dlMP x) 


where A and B are dehned in (21) and (23) respectively. 


The constants A and B can be simplihed to: 


A = 


[APA 

l-oc (^2{y) 


(^1 - (^i'^)'x{dy) dy, 


and 


B = 


(p + hVl -P^) 


Ky) 


c^iiy)(^2iy) 


{ai - ai^)n{dy) dy. 


A- 


do 


On substituting the derivatives of Pq in (|40|), we get 

/e = Ui + y/e 

= ai \fk 


A 


B 

+ ~ 
CTl. 

X 

Wr 


B-A/2 


+ o(v^) 


(39) 


(40) 


(41) 


(42) 


(43) 


We see that, as in the iio-arbitrage pricing case (see 0), the corrected implied volatility is 
an affine function of log-moneyness-to-maturity ratio (LMMR), where 


LMMR = 


\n{strike price/asset price) —x 


as 


time to maturity t 

For calibration purposes, it is convenient to rewrite the formula for implied volatility, R, 

h = a{LMMR) + d + o(^/e), (44) 
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where a = ^ 3 ^ , 


and d = ai + ^ (^B — . On calibration of a and d from market 

implied volatility data, we can determine A and B from the formulas A = and 


0-1 


In [1] an affine function of LMMR is fitted to S&P 500 European call option implied 
volatility data and a and d are estimated to be -0.154 and 0.149 respectively. Using these 
estimates we graph the implied volatility surface given by (44), see Figure 

Taking Yf to be an OU-process and (Ji{y) an arctangent volatility function, we compute 
the implied volatility from the corrected asymptotic option price formula in Theorem 1. To 
be precise, we take = 0.3 — arctan(?/), < 72 ( 1 /) = 0.2, m = 0, b{y) = 1 and p = —0.2 

in our stochastic volatility model ([^. We fix r = 0.25 and assume e = 0.004. Figure j^gives 
the implied volatility as a function of log moneyness i.e. —x = \ogK/S, for three different 
values of the risk parameter r]. The skew of the implied volatility function is clearly seen. 
Observe that as the value of rj increases, the implied volatility gets shifted down. This was 
also observed in [12] (see Figure 3 in |I2])- In Figure implied volatility appears to be an 
almost linearly decreasing function of log moneyness. This agrees with the formally derived 
formula for corrected asymptotic implied volatility in 


Remark 4.1. The corrected implied volatility formula only depends on the volatility risk pa¬ 
rameter 7 ] and the correlation term p, and has no dependence on the risk aversion parameter 
7. By calibration, the correlation coefficient p and the risk parameter p can he obtained from 


(41) and (42) 
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Figure 1: Implied volatility surface given by (44) with a 


-0.154 and d = 0.149. 



Figure 2: Implied volatility as a function of log moneyness for the arctangent stochastic 
volatility model, when r] = —0.25, 0, 0.25. 
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A Appendix 

Lemma A.l. Suppose f G C'ft(]R) is a bounded continuous function which is centered with 

respect to the invariant distribution tt, i.e. J fdn = 0. Then, the Poisson equation 


Bv = f, 


has a solution that grows at most logarithmically and has bounded derivative v'. 

Proof. We construct a solution that satisfies the required growth condition. 

Suppose / is bounded and n is a solution of Bv = f. Recall the definition of the 
differential operator H in ([^. Multiplying the equation, Bv = /, by the integrating factor 
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Without loss of generality, we can assume m = 0, so 




v'(y) = JJJe° ^2 

(^Jy) 


Integrating, we get 


-dz 


J-00 (^Ju) 


(A.l) 


ffj'^f^dz 

under the assumption that hmy^_oo e ‘" 2 (^) v [y) is 0. Let us denote the right hand side of 


(A.l) by the function G{y) i.e. 

G{y) := e^° '^2^ 


^0 J JfJ) Jo Jjj, 


Loo crU) 


e du. 


Using I’hopitars rule as ||/| —>■ 00 and the boundedness of / and ct 2, we see that G is a 
bounded function. Then viy) := ^ G{u)du gives us a solution to the Poisson equation 
Bv = f. We will show that this solution v to the Poisson equation has at most logarithmic 
growth. 

Recall the invariant measure tt given in ([^. Using the centering condition J /dvr, we get 

v^(y) = r 

Jy 

Let y > 1. By the boundedness of / and (T2, we can bound 


rv 2z j_ 

Jo 


|^^'(|/)| <ce° LW 


2 /o 


dz 


ly 


e du 


(for some constant c > 0) 

Jo f J^ . 


= ce o'- 


-u aju) 


e ""at") du, 


which on integrating by parts gives 


ri 

Jo 


= ce 2 


<J 




1 -fo^G^dz 


-e 

Ly 




1 fo 


-dz 


—e du 

u^ 


y 


(A.2) 


for some Ci > 0. 
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Let y < —1. We repeat the same argument using the bound 


|u'(?/)| < e° 

ry 2z 


■^0 2 |/( m )| ^/ o “ 


ry 

Jo 


= ce 
C2 




dz 


-00 

1 -S^^dz 

-y 


e du 


^ -Jo f 1 

—e ‘^ 2 ^ _ / e 0 <^ 2 W du 

1 -^ 


(A.3) 


< 


-y 


for some C 2 > 0. The boundedness of \v'{y)\ for all y together with the bounds (A. 2 ) and 
(A.3) when \y\ > 1 gives us 


|^^(l/)| < C'ilog(l +I 2 /I)+ C 2 , V?/e M, 


for some Ci, C 2 > 0. □ 

Lemma A.2. = 0(—loge). 

Proof. Recall that is chosen to be a positive number such that 

> \u 2 iT,x,y) \ + a/c \u 3 {T,x,y) \ - e{y - mf, 

for all r, x, y. Since U 2 and M 3 are bounded in r and x and have at most logarithmic growth 
in ?/, we can write 


\u 2 {T,x,y) \ + a/c \u^{T,x,y)\ < ci log(l + {y - mf) + C 2 for all T,x,y, 

for some positive constants Ci and C 2 . So it suffices to choose such that 

> Cl log(l + {y- mf) + C 2 - e{y - mf, (A.4) 

for all T,x,y. The quadratic term {y — m)^ grows faster than log(l + {y — m)f for large 
\y\, therefore the maximum of the right hand side in (A.4) is attained at finite values of y. 
We will determine the maximum value of the right hand side of (A.4) and choose to be 
larger than that. 

Let ?/o denote a point of maximum of the right hand side of (A.4), then 
^ [ci log(l + {y- mf) + C 2 - e{y - m)f \y=y^ = 0 , 


i.e. 


ci 2 (|/o - m) 

1 + ( 1/0 - mf 


e2{yo - m) = 0 . 
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Solving for we get, either j/o 


and the right hand side of (A.4) 


Therefore, for small enough e, it 


: m, in which case choose > C 2 or 
{yo - m)2 = y - 1, 
becomes 

Cl log(y) + C 2 - Cl + e. 
suffices to choose = — 2ci loge. 


□ 
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